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Resume Os problemas de optimizagao dinamica (em espagos de fungoes) trata- 
dos pelo calculo das variagoes, sao normalmente resolvidos por recurso as condigoes 
necessarias de Euler-Lagrange, que sao equagoes diferenciais de segunda ordem (ou 
de ordem superior, quando os problemas variacionais envolvem derivadas de ordem 
superior a um). Estas equagoes sao, em geral, nao lineares e de diflcil resolugao. 
Uma forma de as simplificar consiste em obter leis de conservagao: primeiros inte- 
grais das equagoes diferenciais de Euler-Lagrange. Os primeiros integrais permitem 
baixar a ordem das equagoes e, em casos extremos, com um numero suficientemente 
grande de primeiros integrais independentes, resolver o problema por completo. Se 
em areas como a Ffsica e a Economia a questao da existencia de leis de conservagao 
e resolvida de forma bastante natural, a propria aplicagao sugerindo as leis de con- 
servagao (e.g. conservagao de energia, conservagao da quantidade de movimento, 
conservagao do rendimento, etc.), de um ponto de vista estritamente matematico, 
dado um problema do calculo de variagoes, o processo de obtengao das leis de con- 
servagao ou, ate mesmo, a demonstragao de que elas existem (ou nao), deixa de ser 
uma questao obvia. Neste trabalho mostramos como um sistema de computagao 
algebrica como o Maple pode ser muito litil na abordagem a estas questoes. Ap- 
resentamos um conjunto de facilidades computacionais simbolicas que permitem, 
de uma forma sistematica e automatica, identificar as leis de conservagao de uma 
dada funcional integral do calculo das variagoes. O algoritmo usado tern por base o 
celebre teorema de Emmy Noether, que associa a existencia de leis de conservagao 
as propriedades de invariancia do problema (a existencia de simetrias variacionais) . 
Varios exemplos ilustrativos sao apresentados, mostrando a utilidade das ferramen- 
tas desenvolvidas. 

Abstract. The dynamic optimization problems treated by the calculus of varia- 
tions are usually solved with the help of the 2nd order Euler-Lagrange differential 

'Partially presented at XXVII CNMAC (Congresso Nacional de Matematica Aplicada e Com- 
putacional), FAMAT/PUCRS, Porto Alegre, RS, Brasil, 13-16 September 2004. Supported by the 
program PRODEP III/5.3/2003 and the R&D unit CEOC. 
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equations. These equations are, generally speaking, nonlinear, and very hard to 
solve. One way to address the problem is to obtain conservation laws of lower or- 
der than those of the corresponding Euler-Lagrange equations. While in Physics 
and Economics the question of existence of conservation laws is treated in a rather 
natural way, because the application itself suggest the conservation laws (e.g., con- 
servation of energy, income/health law), from a strictly mathematical point of view, 
given a problem of the calculus of variations, it is not obvious how one might derive 
a conservation law or, for that matter, if it even has a conservation law. The ques- 
tion we address is thus to develop computational facilities, based on a systematic 
method, which permits to identify functionals that have conservation laws. The 
central result we use is the celebrated Noether's theorem. This theorem links con- 
servation laws with the invariance properties of the problem (with symmetries), and 
provides an algorithm for finding conservation laws. Thus the problem is reduced 
to the one of finding the variational symmetries. We show how a Computer Algebra 
System can help to find the symmetries and the conservation laws in the calculus 
of variations. Several illustrative examples are given. 

Mathematics Subject Classification 2000: 49-04, 49K05, 49S05. 

1. Intro dugao 

O calculo das variagoes e uma area classica da matematica, com mais de tres seculos 
de idade, extremamente activa no seculo XXI, e com inumeras aplicacoes praticas 
na mecanica, economia, ciencias dos materials, ciencias do espago e engenharia. O 
calculo das variagoes esta na origem de muitas areas mais recentes, como sejam a 
analise funcional e o controlo optimo 0J0|S]. 

A computagao algebrica, tambem chamada de computagao cientifica ou simboli- 
ca, permite trabalhar com expressoes matematicas de maneira simbolica, nao numc- 
rica, e e uma area de investigagao moderna, que surgiu na segunda metade do seculo 
XX. Se e verdade que os actuais sistemas de computagao algebrica sao extremamente 
poderosos, tambem nao e menos verdade que muito existe por fazer no sentido 
de se colocarem os computadores a realizar tarefas matematicas verdadeiramente 
interessantes. 

O uso da teoria do calculo das variagoes, na resolugao de problemas concretos, 
requer calculos nao numericos: no calculo das variagoes a presenga de fcnomenos 
como o de Lavrentiev faz com que as solugoes algebricas exactas sejam muito 
mais convenientes do que as solugoes numericas. Tendo em conta que a resolugao 
dos problemas do calculo das variagoes passa quase obrigatoriamente pela resolugao 
das equagoes diferenciais de Euler-Lagrange, equagoes estas, em geral, de dificil 
resolugao, neste trabalho propomos um conjunto de procedimentos computacionais 
algebricos que permitem automatizar o processo de obtengao de leis de conservagao 
(primeiros integrals das equagoes de Euler-Lagrange). Como e bem conhecido da 
teoria das equagoes diferenciais, estas leis de conservagao sao de extrema utilidade, 
permitindo baixar a ordem das equagoes |12| . 

O presente trabalho encontra-se organizado em sete secgoes. Em fazemos 
uma breve apresentagao do calculo das variagoes, da condigao necessaria de Euler- 
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Lagrange (definida em Maple na SjfjQ) e definimos lei de conservagdo. Em il3~l defin- 
imos simetria (invaridncia de urn problema do calculo das variagoes) e apresenta- 
mos um metodo construtivo para as obter (implementado em Maple na por 
intermedio do procedimento Simetria). Em H4.I formulamos o celebre princi'pio de 
Noether, que nos da uma formula explfcita para as leis de conservagao do calculo das 
variagoes, em fungao dos geradores que definem a simetria. O teorema de Noether 
e definido em Maple na por intermedio do procedimento Noether. Varios ex- 
cmplos, mostrando o modo de uso e a utilidade das fungoes Maple definidas em 
SOO s &o apresentados em Por fim, apresentamos em 371 aleumas conclusoes e 
perspectivas de trabalho futuro. 



2. Calculo das Variagoes 

Neste trabalho consideramos problemas do calculo das variagoes de ordem superior: 
minimizar uma funcional integral 

J[x(-)]= / L(t,x(t),x(i),...,xM( t )) dt) (2.1) 



onde o Lagrangeano L e uma fungao real que assumimos ser continuamente difer- 
enciavel em [a, b] x ]^ nx ( m + 1 )- t <eR, e & variavel independente; x(t) = [xi(t) X2{t) 
■ ■■x n (t)] T e R™, as variaveis dependentes; x«(*) = [%& . . . %^] T e 

R™, com i = 1, . . . , m, as derivadas de ordem i das variaveis dependentes em ordem 
at; e x(t) = x^'(t). Assim, o Lagrangeano considerado, depende explicitamente 
de uma variavel independente, de n variaveis dependentes e das suas m primeiras 
derivadas. Para m = 1 obtemos o problema fundamental do calculo das variagoes. 

Na minimizagao da funcional 1|2.1|) e usual recorrer-se ao sistema de equagoes de 
Euler-Lagrange 

dL _ d l / dL 



=i 



onde -Mir = \-^iTT -^77T • ■ • -^tttI e L, e suas derivadas, sao avaliadas ao longo de 
(i,x(t),x(f),...,xH(f)). 

Definigao 2.1. v4s solugoes da equagdo de Euler-Lagrange chamamos ex- 

tremals. 

Observagao 2.1. Para o 'problema fundamental do calculo das variagoes, i.e., 
quando a funcional nao depende de derivadas de x(t) de maior ordem do que a 
primeira (m — l),o sistema de equagoes de Euler-Lagrange \2.2\i reduz-se a 



dL_d_ (dL 
dx dt \ dx 



0. 



Observagao 2.2. Se a funcional nao envolver mais do que uma variavel dependente 
(n = 1), todas as grandezas presentes em \2.1)) e \2.2)) sao escalares. 
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Em definimos o procedimento EulerLagrange que tern por entrada o La- 
grangeano e como sai'da o sistema de equagoes de Euler-Lagrange correspondente, 
que resulta da aplicagao da equagao l|2.2|) . 

Deflnigao 2.2. Uma fungao t — > </>(i, x(i), x(t), . . . , x( fc ) (t)), fc < 2m, que se man- 
tenha constante ao longo de todas as extremals do problema \2.1\) . 6 chamada de 
primeiro integral (de ordem k) da equagao de Euler-Lagrange A equagao 

^(t,x(i),x(i), . . . ,x (fc) (f)) = const 

chamamos lei de conservagao (de ordem k). 

As leis de conservagao sao muito liteis, pois permitem reduzir a ordem das 
equagoes differencials de Euler-Lagrange. Com um numero suficientemente elevado 
de primeiros integrals independentes, e mesmo posslvel determinar explicitamente 
as extremais. Consideremos, a ti'tulo de exemplo, o seguinte problema (n — m — 1): 

,6 

J[x(-)] = / (x 2 {t) - x 2 {t)) dt — > min . 

J a 

Neste caso o Lagrangeano e dado por L(x, x) = — x 2 + x 2 e a equagao de Euler- 
Lagrange H2.2f> reduz-se a x(t) + x(t) = 0. Se multiplicarmos esta equagao diferencial 
por cos(£) obtemos ^ (icos(i) +a;sin(t)) = 0; enquanto que se a multiplicarmos 
por — sin(t) obtemos 4? (— isin(i) + xcos(t)) = 0. Temos entao duas leis de conser- 
vagao: 

\xcos(t) + xsin(i) = ci , _ 2 ^ 

1 — i sin(t) -I- x cos(t) = C2 . 

Resulta de imediato das duas leis de conservagao (|2.3|l que as extremais tern a forma 
x{t) = ci sin(i) + C2 cos(t). Claro que este exemplo e trivial: a equagao diferencial 
de Euler-Lagrange podc facilmcntc scr resolvida sem recurso a leis de conservagao, 
uma vez que ela e linear. Usando o Maple e o nosso procedimento EulerLagrange 
farfamos: 

> L := v~2 - x~2: 

> dsolve (EulerLagrange (L,t ,x,v) ) ; 

{x (t) = .CI sin (t) + -C2 cos (t)} 

No caso nao linear o comando Maple dsolve nem sempre e capaz de obter solugoes 
explfcitas (e.g. para a equagao nao linear de Emden-Fowler considerada no Exemp- 
lo EH ou para o problema de Kepler tratado no Exemplo 15. 2(1 . Nessas situagoes as 
leis de conservagao sao muito liteis. Coloca-se entao a seguinte questao: Dada uma 
funcional integral do tipo \2.1\l . como obter leis de conservagao? Esta questao foi 
resolvida por Emmy Noether em 1918 |llj : se a funcional for invariante sob deter- 
minado tipo de transformagoes (transformagoes de invariancia ou simetrias) , entao 
existem formulas explicitas para as leis de conservagao. O Teorema de Noether en- 
contra muitas aplicagoes em campos concretos da Engenharia, como a Sismologia e 



Computagao Algebrica no Calculo das Variagdes 



5 



a Metalurgia A dificuldade na sua aplicagao reside na obtengao das simetrias (na 
obtengao das transformagoes de invariancia). Neste trabalho automatizamos, por 
recurso ao sistema de computagao algebrica Maple PP , a determinagao de simetrias 
e a correspondente aplicagao do Teorema de Noether. 



3. Simetrias Variacionais 

Para o estudo das propriedades de invariancia das funcionais do calculo das variagoes 
considera-se uma familia uni-parametrica de transformagoes h s it, x) que forma um 
grupo local de Lie [101 12]. A familia uni-parametrica h s (t, x) representa um conjunto 
de n + 1 transformagoes de [a,b] x M™ em M. x K™ 

t s = hf(t,x), xf =/i*.(i,x), com % = 1, ...,n, (3.1) 

a que correspondem n + 1 geradores infinitesimals representados por 



T(t,x)= |-fc?(t,x) 
as 



s=0 ° S 



com i = 1, . . . ,n. (3.2) 



s=0 



Definigao 3.1. A funcional )2.1)) diz-se invariante no intervalo [a,b] sob as trans- 
formagoes uni-parametricas \3. 1)) se, para todo o s suficientemente pequeno, 

f L(t, x(t), x(i), . . . , x< m > (t)) dt = f L(t s ,x s (t s ), x s (0, x fl < m > (t s j) dt s , 

em qualquer subintervalo [a,f3] C [a, b]; com a s — /i£(a,x(a)) e [3 s = hi (/?, x(/3)). 

Observagao 3.1. iVas condigoes da definicao \3.1\ as transformagoes uni-parame- 
tricas \3.1\l constituem uma simetria variacional da funcional )2.1)) . 

O teorema que se segue estabelece uma condigao necessaria e suficiente de in- 
variancia, de extrema importancia para os objectivos a que nos propomos. 

Teorema 3.1 ([15j). A funcional )2.1)) e invariante sob as transformagoes uni- 
parametricas 13.1)) . com geradores infinitesimais T e X 13.2)) . se, e apenas se, 

i=0 



onde 



p o =Xj p ^ = ^- x (^)^, i = 0,...,m-l. (3.4) 



Em 1)3.3)) e \3.4\j assumimos que T e X = \X\ Xi ■ ■ ■ X n ] sao avaliadas em fungao 
de (t, x) e p l em fungao de (t, x, x, . . . , x^), com i = 0, 1, . . . , m. 
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Corolario 3.1.1. Quando o Lagrangeano L ndo depende de derivadas de x(t) 
ordem superior a primeira (m — 1), a equagdo 13. 3\) toma a forma 



dL 
~dt 



T 



dL 
dx 



dL 
dx 



dX . dT\ 



L- 



dT 
~dt 



0, 



com 



ide 



dT 
~dt 



dT 
~dt 



dT 
dx 



dX 
dx 



dX dX 

dxi dx2 



dX 
~dt 



dX 

dx n 



&X 

dt 



dx 
dx" 



dX 2 0X 2 
dxi 8x2 



dX n dX n 
dx± 0x2 



dXx 

dx. n 



dX„ 

dx u 



Observagao 3.2. Todas as derivadas totais presentes em I^S.S}) e \3.4\j podem ser 
expressas por derivadas parciais, usando as igualdades ^ = ^ + |^ • x e 



dp' 1 
~~dt 



dp 1 
~dt 



dx( k ) 



o,. 



onde 



dp 1 

<9x( fc ) 



dp 1 dp 1 



dx^ dx. 



(fe) 



dp 1 



dx 



(fe) 



dp\ 


dp\ 


dp\ 


dx\ 
dpi 


dpi 


dxi k) 
dpi 


dx[ k) 


Ox^ 


dx^ 


dp z n 






dxf^ 


94 fc) 


dx^ 



O teorema l3.ll para alem de servir de teste a existencia de simetrias, estabelece 
urn algoritmo para a determinagao dos correspondentes geradores infinitesimals. 
Como veremos, este facto e crucial: o teorema de Noether fTeorema l4.1f) afirma que 
as leis de conservagao associadas a uma dada simetria variacional apenas dependem 
dos geradores infinitesimals. 

Dado um Lagrangeano L, determinamos os geradores infinitesimals T e X de 
uma famflia uni-parametrica de transformagoes simetricas pelo seguinte metodo. A 
equagao <|3.3|l e uma equagao diferencial nas n + 1 funcoes incognitas T, X\, X2, ■ ■ . , 
e X n , que pretendemos determinar. Porem, a equagao tern de permanecer valida 
para todos os Xi, i — 1, . . . ,n. Como as fungoes T, X\, X2, . . . , e X n dependem de 
n, ao substituirmos, na equagao (|3.3|) . L c todas as suas derivadas 
parciais pelos seus valores, obtemos um polinomio nas n x m variaveis x\, . . . ,x n , 

(2) (2) (m) (m) ,„ . „ - ,,. , 

x\ , . . . ,Xn , ■ ■ -,Xi ,...,Xn e suas potencias. Para que a equagao sej a valida 
para todos os valores das variaveis do polinomio, todos os seus coeficientes devem 
ser nulos. Notamos que os termos do polinomio poderao ser em maior niimero que 
as incognitas do problema (n + 1), pelo que a condigao necessaria e suficiente (|3.3[) 
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pode conduzir a um sistema de equagoes sem solugao. Tal facto significa apenas 
que nem todas as funcionais integrals do calculo das variagoes admitem simetrias 
variacionais (ver Exemplo l5.5|l . O sistema de equagoes a resolver, para a obtengao 
dos ger adores, e um sistema de equagoes diferenciais as derivadas parciais. No 
entanto, ao contrario das equagoes diferenciais ordinarias de Euler-Lagrange, em 
geral nao lineares e de diffcil resolugao, este sistema e linear em f^-, e 

A resolugao do sistema de equagoes diferencias parciais que deriva da expressao 
(|3.3Jl . nomeadamente quando lidamos com valores denem superiores a unidade, 
envolve um mimero muito elevado de calculos, o que torna premente a necessi- 
dade de nos munirmos de ferramentas computacionais que automatizem o trabalho. 
Com esse fim, desenvolvemos um procedimento em Maple, dcsignado Simetria. O 
respectivo codigo, assim como o dos restantes procedimentos por nos desenvolvidos, 
e apresentado na secgao||0 O procedimento Simetria tern por entrada a expressao 
que caracteriza o Lagrangeano e como sai'da os respectivos geradores infinitesimals. 
No caso do Lagrangeano nao admitir qualquer simetria, obtemos geradores nulos 
(cf. Exemplo l5.5|) . 

Na secgao que se segue mostramos como os geradores, obtidos por intermedio 
do nosso procedimento Simetria, podem ser usados na obtengao explfcita de leis de 
conservagao. 



4. Leis de Conservagao 

Emmy Noether foi a primeira a estabelecer a relagao entre a existencia de simetrias 
e a existencia de leis de conservagao. Esta ligagao constitui um princfpio universal, 
passfvel de ser formulado na forma de teorema nos mais diversos contextos e sob as 
mais variadas hipoteses P 131 IH1 IT31 IHl IT5] . 

Teorema 4.1 (Teorema de Noether [15j). Se a funcional $2.1)) e invariante 
sob as trans formagoes uni-parametricas $S.1\) . com geradores infinitesimals T e X, 
entao 

rn / m \ 

2 ** • p*" 1 + I I, - • X W T = const, te[a,b], (4.1) 

i=l V i=l / 



com 



8L 



Tl _! dL d* 1 . 

d* 1 d& 



E ( x (, + D) 



dt dt ^ V / <9x« 

fe=0 
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onde 

e onde se considera a grandeza ~$> l avaliada em (t, x(£), x(i), . . . , x ( 2m ~ l ) (t)), j = 
l,...,m. 

Corolario 4.1.2. Quando o Lagrangeano L nao depende das derivadas de x(i) de 
maior ordem do que a primeira (m = 1), a equagao 14 reduz-se a 

-— ■ X + ( L — — • x ) T = const. 
ax \ ox J 

As leis de conservagao que procuramos sao obtidas substituindo nas equagSes 
(|3.4|l e (|4.1I) os geradores infinitesimais T e X encontrados pelo metodo descrito na 
secgao anterior. Na seccao HP definimos o procedimento Noether. Estc proccdimcnto 
tem por entradas o Lagrangeano e os geradores infinitesimais, que sao obtidos por 
intermedio do nosso procedimento Simetria, e como saida a correspondente lei de 
conservagao (|4. 1|> . Resumindo: dado um problema do calculo das variagoes (|2. Ill . 
obtemos as leis de conservagao, de uma forma automatica, atraves de um processo 
de duas etapas: com o nosso procedimento Simetria obtemos todas as possiveis 
simetrias do problema; recorrendo depois ao nosso procedimento Noether, baseado 
no teorema |4.1l obtemos as correspondentes leis de conservagao. Na secgao seguinte 
apresentamos alguns exemplos que ilustram todo o processo. 

5. Exemplos Ilustrativos 

Consideramos agora varias situagoes concretas, mostrando a funcionalidade e a 
utilidade das ferramentas desenvolvidas. 

Exemplo 5.1. Comegamos com um exemplo muito simples em que o Lagrangeano 
depende apenas de uma variavcl depcndcntc (n = 1) c nao existem derivadas de 
ordem superior a primeira (m = 1): L(t,x,x) = tx 2 . 
Com a definigao Maple 

> L:= t*v~2; 

L := tv 2 

o nosso procedimento Simetria detcrmina os geradores infinitesimais das simetrias 
do problema do calculo das variagoes em consideragao: 

> Simetria(L,t,x,v) ; 



dx( fe ) 



a** 

dx[ k > 

a** 
a4 fc) 




a** 
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{T (t, x) = (2 .CI In (t) + _C3) t, X (t, a:) = _Ci x + .C2} 

A fanrilia de geradores depende de tres parametros que advem das constantes de 
integragao. Com as substituicoes 1 

> subs (_C1=1 , _C2=0 , _C3=0 ,'/„); 

{T (t, x) = 2 In (t) t, X (t, x) = x} 

obtemos os geradores descritos em pp. 210 e 214]. A lei de conservagao correspon- 
dente a estes geradores e facilmente obtida por intermedio do nosso procedimento 
Noether: 



> LC := Noether (L, t ,x, V,*/.) : 



d ,.s ,2 ( d 
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LC := x it) t—x (t) - t l—x {t)j In (t) = const 

E, neste caso, muito facil verificar a validade da lei de conservagao obtida. Por 
definigao, basta mostrar que a igualdade e verificada ao longo das extremais. A 
equagao de Euler-Lagrange e a equagao diferencial de 2 a ordem 

> EulerLagr ange (L , t , x , v) ; 

-*!«(*)-24*(t)=0 
e as extremais sao as suas solugoes: 

> dsolve ('/.); 

{x(t) = -CI +_G21n(f)} 

Substituindo as extremais na lei de conservagao, obtemos, como esperado, uma 
proposigao verdadeira: 

> expand(subs (°/,,LC) ) ; 

-C2.C1 = const 

Exemplo 5.2 (Problema de Kepler). Analisamos agora as simetrias e leis de 
conservagao do problema de Kepler [21 p. 217]. Neste caso o Lagrangeano tern duas 
variaveis dependentes (n — 2) e nao envolve derivadas de ordem superior (m = 1): 

L(t,q,q) = -(^ + gi) 1 



2 v ' 

Vamos determinar a formula geral das leis de conservagao. Neste caso nao e possivel 
validar a lei de conservagao por aplicagao directa da definigao, como fizemos para 
o exemplo anterior, pois o Maple nao e capaz de resolver o respectivo sistema de 
equagoes de Euler-Lagrange 

> L:=m/2*(v[l] ~2+v[2] ~2)+K/sqrt(q[l] ~2+q[2] ~2) ; 



1 sinal de percentagem (%) e um operador usado em Maple para referenciar o resultado do 
comando anterior. Para uma introdugao ao Maple, veja-se 0. 
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1/2 m(V +v 2 2 ) + K 



> EulerLagrange(L,t, [q[l] ,q[2]] , [v[l] ,v[2]] ) ; 

f d 2 ,s Kqi (t) „ d 2 , . ATa 2 (t) 

-m — 91 (t) - ,2,3/2 - °' "™^2 (t) - ' 2 ,3/2 = 

I d * (?i (*) + 92 (t) 2 ) di (gi (t) 2 + q 2 (t) 2 ) 

> Simetria(L, t, [q[l],q[2]], [v [1] , v [2] ] ) ; 

{X 1 (t, qi ,q 2 ) = -C2q 2 , T(t,qi,q 2 ) = -Cl, X 2 (t, qi, q 2 ) = --C2 qi} 

> Noether(L, t, [q[l],q[2]], [v[l],v[2]], "/.) : 

> expand (°/ ) ; 

7 7 / 7 \ 2 

.C2q 2 {t)m—q l {t)--C2q 1 {t)m—q 2 {t)- 1/2 .CI f^?i(*)J ™ 

, /„ ™ ( d , X 1 -C1K 
-1/2 _ CI — q 2 ft) mH = = const 

J ^ qi{t f + q2 (tf 

Exemplo 5.3. Vejamos o caso de um Lagrangeano com duas variaveis dependentes 
(n = 2) e com derivadas de ordem superior (m = 2): 

L(t, x, x, x) — x\ + x 2 

> L:=v[l] ~2+a[2] ~2; 

T 2.2 

L :— vi + a 2 

> SimetriaCL, t, [x[l],x[2]], [v[l],v[2]], [a[l] ,a[2] ] ) ; 

{ T(t,xi 7 x 2 ) =-Clt + _C2, 

Xi (t,xi,x 2 ) = 1/2 .CI X-L+.C5, X 2 (t,xi,x 2 ) = S/2.C1 x 2 + .C3t + .C4 } 

> LC := Noether(L, t, [x[l],x[2]], [v[l],v[2]], [a[l],a[2]], "/,) ; 



LC ~ 2 (1/2 .CI Xl ft) + _C5) —xi (t) - 2 (3/2 _Ci x 2 (t) + .C3 1 + _C4) —x 2 (t) 

at at J 

+ 2 (_C9 + 1/2 

d , \ 2 ( d 2 , A 2 „ d , , d 3 



dt' 1 ®) -^WJ +2 5 « a (*) sr « a (t)j(_<7i* + _G») = »™t 

Tal como para o Exemplo 15.11 tambem aqui e facil verificar, por aplicagao directa 
da definigao, a validade da lei de conservagao obtida: 

> EulerLagrangeCL.t, [x[l] ,x [2] ] , [v [1] , v [2] ] , [a[l] ,a[2]] ) ; 
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> dsolveC/.) ; 

{ xi (t) = .Clt + .C2, x 2 it) = 1/6 .C3t 3 + 1/2 _C4 t 2 + .C5t + .C6 } 

> expand(subs (°/ t> LC) ) ; 

2 .01 .05 - 3 .C3 .CI .06 + .CI .05 .04 - .C4 2 .C2 + 2 _C3 .C5 .02 = const 

Exemplo 5.4 (Emden- Fowler). Consideremos o problema variacional definido 
pelo Lagrangeano 

> L:= t~2/2*(V2-(l/3)*x~6) ; 

L - 2 

A respectiva equagao diferencial de Euler-Lagrange e conhecida na astrofisica como 
a equagao de Emden- Fowler p. 220]: 

> EL := EulerLagrange(L,t,x,v) ; 



2 



EL := {-2*|U (t) - t 2 ^x (t) - t 2 (x (t)) 5 = 

Encontramos os geradores infinitesimals, que conduzem a uma simetria variacional 
para a funcional de Emden-Fowler, por intermedio da nossa fungao Simetria: 

> S := Simetria(L,t,x,v) ; 

{t C1 
X(t,x) = ±-,T(t,x) = _Clt 

Por exemplo, 

> S2 := subs(_Cl=-6,S) ; 

52 := {T(t, x) = -&t,X(t, x) = 3x} 

Aplicando o Teorema de Noether (Teorema I4.1(l . estabelecemos a seguinte lei de 
conservagao: 

> simplify (Noether (L,t ,x,v, S2) ) ; 



const 

Exemplo 5.5 (Thomas- Fermi). Mostramos agora um exemplo de um problema 
do calculo das variagoes que nao possui nenhuma simetria variacional. Seja 

> L:=l/2 * v"2 + 2/5 * (x~ (5/2) )/ (sqrt (t) ) ; 

v 2 2x2 



L~— + 



2 5Vt 

A equagao de Euler-Lagrange associada a este Lagrangeano corresponde a equagao 
diferencial de Thomas-Fermi [2| p. 220]: 
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> EL := EulerLagrange(L,t,x,v) ; 



EL := 




Vt 



= o 



A nossa fungao Simetria devolve, neste caso, geradores nulos. Como explicado na 
E I3.I isto significa que este problema do calculo das variagoes nao admite simetrias. 

> Simetria(L, t, x, v) ; 



Exemplo 5.6 (Oscilador Harmonico com Amortecimento). Consideremos 
um oscilador harmonico com forga de restituicao — kx, submerso num lfquido de tal 
modo que o movimento da massa m e amortecido por uma forga proportional a sua 
velocidade. Recorrendo a segunda lei de Newton obtem-se, como equagao de movi- 
mento, a equacao diferencial de Euler-Lagrange associada ao seguinte Lagrangeano 
EH pp. 432-434]: 

> L:=l/2 * (m*v~2-k*x~2)*exp((a/m)*t) ; 



' 2 v ; 

Para determinar um primeiro integral da equacao de Euler-Lagrange, encontramos 
os geradores T e X sob os quais a funcional integral J[x(-)] = J Ldt e invariante: 

> Simetria(L, t, x, v) ; 



{X(t,x) = Q,T(t,x) = 0} 



A fungao Noether resulta num trufsmo: 



> Noether (L, t, x, v, '/„) ; 



= const 





> S:= subs(_Cl=l,7.) ; 



S:= {x(t,x) = -^,T(t,x) = l} 




> simplify (Noether (L, t, x, v, S)); 
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6. Procedimentos Maple 

Os procedimentos Simetria, Noether e EulerLagrange, descritos e ilustrados nas 
secgoes anteriores, sao agora definidos usando o sistema de computagao algebrica 
Maple (versao 9). 

Simetria determina as simetrias de urn Lagrangeano de varias variaveis depen- 
dentes e com derivadas de ordem superior, de acordo com a seccao I3~l 

Devolve: 

- conjunto/vector de geradores infinitesimals das transformagoes simetricas. 
Forma de invocagao: 

- Simetria(L, t, x, xl, x2, xm) 

Parametros: 

L - expressao do Lagrangeano; 
t - nome da variavel independente; 

x - nome, lista de nomes ou vector de nomes das variaveis dependentes; 

xi - nome, lista de nomes ou vector de nomes das derivadas de ordem i das 

variaveis dependentes; 

Simetria: =proc(L : : algebraic , t : :name , xO : : {name , list (name) , 'Vector [column] ' 

(name) } ,xl : : {name , list (name) , 'Vector [column] ' (name) }) 
local n,m,xx,P ,EqD,SysEqD, Sol,xi ,Tdt , soma, V,r , i , j ,k; 
if nargs<4 then print('N° de args insuf iciente . ' ) ; return; 
elif not type ( [args [3 .. -1] ],{' list ' (name) ,' listlist ' (name) , 

'list ' ( 'Vector [column] ' (name) ) }) 
then print ('Erro na lista das var. depend, ou suas derivadas.'); return; 
end if ; 

unassign( 'T' ) ; unassign( 'X' ) ; 

xx :=convert (xO, ' list ' ) [] ; n:=nops( [xx] ) ; m:=nargs-3; 

xi : = [seq( Vector (convert (args [i] , 'list')) , i=3 . .m+3)] ; 

Tdt :=diff (T(t ,xx) ,t)+Vector [row] ( [seq(diff (T(t,xx) ,i) ,i=xx)] ) .xi [2] ; 

if n>l then P := [Vector ( [seq(X [i] (t ,xx) , i=l . .n)] )] ; 

else P:=[Vector([X(t,xx)])] ; end if; 

for i from 1 to m do 

V:=Vector(n) : 

for k from to i-1 do 

V:=V+Matrix(n, (r , j ) ->dif f (P [i] [r] , xi[k+l] [j] )) .xi[k+2] ; 

end do : 

P: = [P[] ,map(diff ,P[i] ,t)+V-xi [i+1] *Tdt] : 
end do : 
soma:=0 : 

for i from to m do 

soma:=soma+Vector [row] ( [seq(dif f (L,r) ,r=convert (xi [i+1] , list) )] ) .P [i+1] : 
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end do : 

EqD : =dif f (L,t)*T(t, xx) +soma+L*Tdt ; 

EqD : =collect (EqD, [seq(convert (xi [i+1] , ' list ' ) [] , i=l . .m)] .distributed) ; 
SysEqD:={coeff s(EqD, [seq(convert (xi [i+1] , 'list') [] ,i=l . .m)])>; 
Sol : =pdsolve (SysEqD, {T(t ,xx) } union convert (P [1] ,' set ')) ; 
if type (xO, 'Vector ' ) then 

return (subs (Sol ,T(t ,xx) ) .Vector (subs (Sol ,P [1] )) ) ; 
else return Sol; end if; 
end proc : 

Noether dados os geradores infinitesimals, determina a Lei de Conservagao de 
urn Lagrangeano de varias variaveis dependentes e com derivadas de ordcm 
superior, de acordo com a seccaoPTl 

Devolve: 

- lei de conservagao. 
Formas de invocagao: 

- Noether(L, t, x, xl, x2, xm, S) 

- Noether(L, t, xs, xls, x2s, xms, T, X) 

Parametros: 

L - expressao do Lagrangeano; 

t - nome da variavel independente; 

x - nome ou lista de nomes das variaveis dependentes; 

xi - nome ou lista de nomes das derivadas de ordem i das variaveis depen- 

dentes; 

S - conjunto de geradores infinitesimals das simetrias (output do procedi- 
mento Simetria); 

xs - vector de nomes das variaveis dependentes; 

xis - vector de nomes das derivadas de ordem i das variaveis dependentes; 

T - gerador da transformagao para a variavel independente (t); 

X - vector com os geradores das transformagoes para as variaveis dependentes 
(xs). 

Noether :=proc (L : : algebraic , t : :name ,x0 : : {name , list (name) , 'Vector [column] ' 

(name) } ,xl : : {name , list (name) , 'Vector [column] ' (name) }) 
local xx,n,m,P ,psi ,LC,xi ,Tdt , V,r , i , j ,k; 

if type (xO, 'Vector ' ) then m:=nargs-5; else m:=nargs-4; end if; 

if m<l then print('N° de args insuf iciente . ' ) ; return; 

elif not type ( [args [3 .. 3+m] ],{' list ' (name) ,' listlist ' (name) , 

' list ' ( 'Vector [column] ' (name) ) }) then 
print ('Erro na lista das var. depend, ou suas derivadas.'); return; 
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elif (type(xO, 'Vector ' ) and not (type (args [-1] , 'Vector [column] ' ) and 

type (args [-2] , algebraic) ) ) or (not type(xO, 'Vector ' ) and 

not type (args [-1] ,' set ') ) then 
print (' Conj . de gerad. invalido.'); return; 
end if ; 

xx :=convert (xO, ' list ' ) [] ; n : =nops ( [xx] ) ; unassign( 'T' ) ; unassign( 'X' ) ; 

xi := [seq( Vector (convert (args [i] , 'list ')) , i=3 . .m+3)] ; 

xi : = [xi [] , seq( Vector ( [seq(x I I i [k] ,k=l . . n)] ) , i=m+l . . 2*m-l)] ; 

Tdt:=diff (T(t,xx) ,t)+Vector [row] ( [seq(dif f (T(t ,xx) ,i) ,i=xx)]) .xi[2] ; 

if n>l then P : = [Vector ( [seq(X [i] (t ,xx) , i=l . .n)] )] : 

else P:=[Vector([X(t,xx)])] : end if: 

for i from 1 to (m-1) do 

V:=Vector(n) : 

for k from to i-1 do 

V:=V+Matrix(n, (r, j)->diff (P [i] [r] ,xi[k+l] [j] )) .xi[k+2] end do: 

P: = [P[] ,map(diff ,P[i] ,t)+V-xi [i+1] *Tdt] : 
end do: 

psi : = [Vector [row] ( [seq(dif f (L, i) , i=convert (xi [m+1] , ' list '))])] : 
for i from m by -1 to 2 do 
V: =Vector [row] (n) : 
for k from to 2*m-i do 

V : =V+LinearAlgebra [Transpose] (xi[k+2]) .Matrix(n, (j,r)-> 

diff (psi[l] [r] ,xi[k+l] [j])) ; end do: 
psi : = [Vector [row] ( [seq(dif f (L, i) , i=convert (xi [i] , ' list ' ) )] ) 

-map(diff ,psi[l] ,t)-V,psi[]] : 

end do : 

LC:=sum('psi[i] .P[i] ' , 'i'=l. .m) + (L-sum( 'psi [i] .xi[i+l] ' , 'i'=l. .m))*T(t,xx) 

=const : 

if type (xO, 'Vector ' ) then 

LC:=eval(LC, [T(t ,xx)=args [-2] ,seq(P[l] [i] =args [-1] [i] ,i=l. .n)]) ; 
else LC : =eval (LC ,args [-1] ) ; end if; 
LC:=subs({map(i->i=i(t) , [xx] ) []>,LC) ; 

LC:=subs({seq(seq(xi[k+l] [i] =dif f (xi [1] [i] (t) ,t$k) ,i=l. .n) ,k=l. .2*m-l)}, 

LC) ; 

return LC; 
end proc : 

EulerLagrange constroi o sistema de equacoes de Euler-Lagrange (|2.2|l , dado um 
Lagrangeano de varias variaveis dependentes e com derivadas de ordem supe- 
rior. 

Devolve: 

- conjunto/vector de equacoes de Euler-Lagrange. 
Forma de invocacao: 

- EulcrLagrange(L, t, x, xl, x2, xm) 
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Paramctros: 

L - expressao do Lagrangeano; 
t - nome da variavel independente; 

x - nome, lista de nomes ou vector de nomes das variaveis dependentes; 

xi - nome, lista de nomes ou vector de nomes das derivadas de ordem i das 

variaveis dependentes; 

EulerLagrange :=proc (L : : algebraic ,t : :name ,x0 : : {name , list (name) , 

'Vector [column] ' (name) } ,xl : : {name , list (name) , 'Vector [column] ' (name)}-) 
local xx,n,m,Lxi,xi,V,EL,i, j ,k; 

if nargs<4 then print('N° de args insuf iciente . ' ) ; return; 
elif not type( [args [3 .. -1] ] ,{'list ' (name) ,' listlist ' (name) , 

' list ' ( 'Vector [column] ' (name) ) }) 
then print ('Erro na lista das var. depend, ou suas derivadas.'); return; 
end if ; 

xx:=convert(xO, 'list') [] ; n : =nops ( [xx] ) ; m:=nargs-3; 
xi:=[seq(Vector (convert (args [i] , 'list')) ,i=3. .m+3)] ; 
V: = [0$n] ; 

for i from 1 to m do 

Lxi: = [seq(diff (L,k) ,k=convert (xi [i+1] , 'list'))] : 
Lxi:=subs({map(k->k=k(t) , [xx] ) []>,Lxi) ; 

Lxi:=subs({seq(seq(xi[k+1] [j] =dif f (xi [1] [j] (t) ,t$k) ,j = l. .n) ,k=l. .m)}, 

Lxi) ; 

V:=V+(-l)~i*map(diff ,Lxi,t$i) ; 
end do: 

EL: = [seq(diff (L,k) ,k=convert (xi [1] , 'list'))] ; 
EL:=subs({map(k->k=k(t) , [xx] ) []},EL) ; 

EL:=subs({seq(seq(xi [k+1] [j] =dif f (xi [1] [j] (t) ,t$k) ,j = l. .n) ,k=l. .m)},EL); 
EL:=EL+V; 

if type (xO, 'Vector ' ) then return convert (map(i->i=0 , EL) , 'Vector [column] ') ; 
elif type (xO ,' list ' ) then return convert (map(i->i=0, EL) ,' set ') ; 
else return op(EL)=0; end if; 
end proc : 

7. Conclusao e Trabalho Futuro 

Os sistemas actuais de computagao algebrica colocam a nossa disposigao ambientes 
de computagao cientifica extremamente sofisticados e poderosos. Disponibilizam 
ja muito conhccimcnto matcmatico e permitem estender esse conhecimento por in- 
tcrmcdio de linguagens de programagao de muito alto nivel, expressivas e intuitivas, 
proximas da linguagem matematica. Tais ambientes permitem a realizagao de uma 
miriade de calculos matematicos simbolicos, com extrema eficiencia, e a definigao 
celere de novas funcionalidades. Neste trabalho usamos o sistema de computagao 
matematica Maple 9 para definir novas fungoes que permitem a determinagao au- 
tomatica de simctrias c lcis de conservagao no calculo das variagoes. Mostramos de- 
pois, por meio de exemplos concretos, como as novas funcionalidades sao de grandc 
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utilidade pratica. Como trabalho futuro pretendemos estender o nosso "package 
Maple" ao caso discreto e aos problemas mais genericos tratados pelo controlo 
optimo. 

Uma abordagem analoga a que fizemos pode ser efectuada para problemas do 
calculo das variagoes discretos no tempo. No caso discreto, o problema fundamental 
do calculo das variagdes consiste em determinar uma sequencia finita x(fc) € R", 
k = M, . . . , M + N, de modo a que a fungao de custo discreta 

M+N-l 

J[x(-)]= ]T i(fc,x(fc),x(fc + l)) 

k=M 

seja minimizada (ou maximizada). Neste tipo de problemas consideramos um in- 
tervalo de tempo de N perfodos, com im'cio num pen'odo fixo M , em que k G Z, 
k = M, . . . ,M + N— 1, e a variavel discreta tempo, e assumimos que o Lagrangeano 
L e continuamente diferenciavel em {M, . . . , M + N — 1} x M™ x R n . A semelhanga 
do que fizemos para o caso continuo, e natural considerarem-se problemas de opti- 
mizagao em que o Lagrangeano envolve diferengas finitas de ordem superior, 

M+N-l 

J[x(-)]= Y. £(fc,x(fc),x(fc + l),...,x(fc + m)), (7-1) 

k=M 

com m > 1, e considerando ainda o Lagrangeno continuamente diferenciavel relati- 
vamente a todos os seus argumentos. Uma condigao necessaria para que x(fc) seja 
uma extremal de (|7.1II e dada pela equagao de Euler-Lagrange discreta 

m dL 

^2g^j(k + m-j,x{k + m-j),...,x(k + 2m-j)) = 0. (7.2) 

3=0 X 

No caso m = 1, a equagao (|7.2|l reduz-se a equagao discreta de Euler-Lagrange 
dL dL 

— (k + 1, x(fc + 1), x(fc + 2)) + g^(k, x(fc), x(fc + 1)) = 0. 

Definigao 7.1. Um problema discreto expresso pela fungao de custo |7. 1)) diz-se 
invariante sob as transformagdes uni-parametricas h s (fc,x), com h°(fc,x) — x, se, 
para todo o s suficientemente pequeno e qualquer que seja o k, 

L(k, x(fc), x(fc+l), . . . , x(fc+m)) = L(k, h s (k, x(fc)), . . . , h s (fc+m, x(fc+m))). (7.3) 

A partir desta definigao conseguimos estabelecer uma condigao necessaria e sufi- 
ciente de invariancia. 

Teorema 7.1. A fungao de custo J 7. J| ) e invariante sob as transformagdes uni- 
parametricas h s (fc,x), com geradores infinitesimals X(fc,x) = ^-h s (fc, x)| s Q , se, 
e apenas se, 

m BT 

]T — •X(fc+ J ,x(fc+ J )) = 0, (7.4) 
j=o x 
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onde o Lagrangeano discreto L e avaliado em (fc, x(fc), x(k + 1), . . . , x(fc + to)). 
Observagao 7.1. Quando m — 1, a equagao j7.4\ ) reduz-se a 

|| • X(fc, x(fc)) + ^ ■ X(fc + 1, x(fc + 1)) = 0. 
Demonstragao (Teorema \7.1\ ). Da condigao de invariancia Ij7.3|l podemos escrever: 

-^-L(k, h s (k, x(fc)), . . . , h s (k + to, x(fc + to))) = . 
ds 



Derivando, obtemos, para s — 0, a igualdade (|7.4|l : 

= 



3=0 



s=0 



com L avaliado cm (k, x(/c), x(fc + 1), . . . , x(fc + m)). 



□ 



E possrvel deduzir um teorema analogo ao teorema de Noether, para o caso discreto, 
que permite estabelecer leis de conservagao discretas (cf. |14|). 

Teorema 7.2 (Teorema de Noether Discreto). Se a fungao de custo \7. l\j 
e invariante, no sentido da definigao \7.1\ sob as transformagoes uni-parametricas 
h s (fc,x) 7 com geradores infinitesimals X(fc,x) = ^-h s (A;,x)| g , entao todas as 
solugdes x(fc) da equagao de Euler- Lagrange \7.<ty satisfazem 

m — 1 

* J (k) • X(fc + j, x(fc +j)) = const (7.5) 

3=0 

onde 

*°(fe) = |^(fe I x(fc),...,x(fc + m)) 
ox 

dL 

= + 1) + (fc,x(fc), . . . ,x(fc + m)), para j = 1, 2, . . . , m - 1 . 

ox-? 

Observagao 7.2. Caso m — 1, a lei de conservagao |7.5| ) reduz-se a 
dL 

-^-{k, x(fc), x(fc + 1)) • X(fc, x(fc)) = const. 

Para problemas discretos no tempo, facilmente se definem em Maple procedimentos 
que implementem as condigoes (17.21 . (|7.4f) e <|7.5ll . 

O controlo optimo pode ser encarado como uma extensao natural do calculo 
das variagoes El ■ A resolugao de problemas do controlo optimo passa normal- 
mente pela aplicagao do Princi'pio do Maximo de Pontryagin, que constitui uma 
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generalizagao das condigoes classicas de Euler-Lagrange e de Weierstrass do calculo 
das variacoes. Em termos praticos, algoritmicos, temos de resolver urn sistema de 
equagoes diferenciais, chamado de sistema Hamiltoniano, depois de eliminados os 
parametros de controlo por intermedio da condigdo de mdximo. Tal como no calculo 
das variagoes, tambem no controlo optimo as equagoes diferenciais ordinarias obti- 
das sao em geral nao-lineares e de dificil resolugao (podem, inclusive, nao ser in- 
tegraveis). As leis de conservagao sao usadas para simplificar essas equagoes e, 
mais uma vez, a questao que se coloca e saber como determinar tais quantidades 
preservadas. Resulta que os resultados classicos de Emmy Noether podem ser gen- 
eralizados ao contexto do controlo optimo, reduzindo o problema ao da descoberta 
de grupos de transformagoes uni-parametricas que deixem o problema de contro- 
lo optimo invariante ^H]- Seria de grande utilidade pratica dispor de facilidadcs 
computacionais que permitissem identificar as simetrias dos problcmas de controlo 
optimo [SJ|2j- 

Serao estas algumas das direcgoes do nosso trabalho futuro. 
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